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Abstract 

In this paper, we introduce a new multivariate discrete distribution which called 
multivariate unification of generalied Apostol type distribution (UGAT). Several prop¬ 
erties are studied as, moments, probability generating function and other properties. 
Also, reliability study of distribution are introduced. Maximum likelihood method is 
used to estimate parameters and numerical method is used to obtain (MLEs). 
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1 Introduction 


Recently, EL-Desouky et al. [5] introduced unification of generalized Apostol Euler, 
Bernoulli and Genocchi polynomials and obtained its explicit formula as 


mTGT;®) 


(_l)r2r(i-fc)„! “ {aoY^{aiY^...{ar-iY- 

{n-rk)i (pTATTTTTT^' 


( 1 . 1 ) 


In the present paper we state a new multivariate discrete distribution (UGAT) distribu¬ 
tion using the explicit formula of the unihcation of generalized Apostol Euler, Bernoulli 
and Genocchi polynomials. Several properties of the new distribution (UGAT) have been 
established as the joint cumulative distribution function, moments, marginal probability 
function and other properties in section 2. Also, in section 2, we present sub models from 
the new distribution (UGAT) 

In 1975, Nakagawa and Osaki [13] were the hrst to study a discrete life time distribution 
which is defined as the discrete counterpart of the usual continuous Weibull distribution. 
In 1982, Salvia and Bollinger [16] introduced basic results about discrete reliability and il¬ 
lustrated them with the simple discrete life distributions. The characterization of discrete 
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distributions has been studied by Roy, Gupta, Gupta in 1999 m- In 1997, Gupta, Gupta 
and Tripahi [7] introduced wide classes of discrete distributions with increasing failure 
rate. In 1997, Nair and Asha |12j introduced some classes of multivariate life distributions 
in discrete time. 

So, In section 3, we introduce some applications of (UGAT) distribution in reliabil¬ 
ity theory. Finally, Section 4 contains the parameter estimation using MLE. In section 5, 
we present a numerical result are obtained using real data. 


2 The model and statistical inference 


In this section, we introduce the multivariate generalized Apostol-Euler Bernoulli and 
Genocchi distribution (UGAT) and we study its structural properties and statistical 
inference. 

Definition 2.1. The discrete random variables Xi, X 2 , ■■■, Xr are said to have UGAT 
distribution with the joint probability mass function 


PyXl — Xi, X 2 — X2j •••? Xj~ — Xr) — / \ 

(n - rk)l{xi +X 2 + ... + Xr + A/3; fe; 

where xi,X 2 , ■■■, Xr = 0, 1, ..., ao,ai, ...a^-i S R-"*") /3 > 0, n € Aq = Au{0} and rk > n. 

( 2 . 1 ) 


We show that many well-known discrete distributions are special cases of UGAT family. 
Some of these are given below 


2.1 Special models 

2.1.1 Lerch distribution 


P{Xi =x) = 


where +l;n + c;p) (£+aP ’ 

£=1 


p 


(x+q^mX (a+l;n+c;p) 


00 


X € A 


2.1.2 

P(Ai 


Hurwitz Lerch zeta distribution 

X__ 

(x+a)^+^Mi^\a+l;n+s+l;0) ’ 

00 


where Mi^\a + T,n + s + l;e) = 

7=1 ^ ^ 


(x G A, s > 0, 0 < a < 1, 0 < 0 < 1), 
see [8]. 
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2.1.3 Good distribution 


P{Xi =x) = 


where n + s + 1; 0) = ^ 


(x)^+iMP(l;n+s+l;0)’ 


(xGA^, 0<6»<1), 


oo 


e 


X 


i=l 


ie) 


s+1 ? 


see 


m- 


2.1.4 Hurwitz Zeta distribution 


PiXi =x) = 


{x+byMi^^ (&;n+CT;l) 


,keN 


where Mi^\b;n + a; 1) = ^ M 


e=o 


2.1.5 Zipf-Mandelbrot distribution 


P{X, = x) =- 


, (x € a > 0, c > 0) 


(x+a)'=M^^^(a+l;n+c;l) 
oo 

where Mi^\a + 1; n + c; 1) = ^ HZ] 


£=1 


2.1.6 discrete Pareto distribution 
1 


P{X, = x) =- 


X € N 


(x)‘=M^^^(l;n+c;l) ’ 
oo 

where m 1^^(1; n + c; 1) = ^ {£+iy ^ HZ] 


e=i 


2.1.7 Geometric distribution 


- 


P{Xl =x) = 


MA^hl;n;p) 

where Mn\l;n;p) = J2 p 


^{1 — p), xeN,{0<p<l) 


£=0 


Lemma 2.1. Let Xi ~ UGAT{ai, (3, k),i = 1, 2,r — 1. Then the marginal joint 
cumulative distribution function 


P{Xi <Xi) = l- (tti-i)- 


+ Xi + l]k] ar) 

A/i’'>(,3;fc;5.) 


( 2 , 2 ) 
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Proof. Since 

P{Xi <Xi) = l- P{Xi > Xi) 


= 1 - 




oo oo 

E E 




{n - rk)\M^n\^] k- a,) (4 + ^2 + ••. + ^* + ... + t " 

£i + l,...£r =0 


Put li — Xi — 1 = I, then 

(l)'-2'-(i-Pn! 


P{X, <Xi) = l- 


E 


(ao)'^(ai)"^...(ai-i)"^+'+h..K-i)'’- 


(n - rk)\Mn^ (/3; k; ar) (^1 + ^2 + ... + ^ + ... + 4 + /? + a^i + 1) 

,h+i,..hr=o 


rk—n 


From fl2.ip . we obtain 


□ 


Theorem 2.1. Suppose Xi,X2, ...,Xr are mutually independent where Xi ~ 
UGAT{ai, /3,k),i = 1,2, — 1. Then the multivariate cumulative distribution 

function is given by 


P{Xi < Xi,X 2 < X 2 , ...,Xr < Xr) = Y\_[^ - 


i=l 


m!{\/ 3 + Xj] k] ar ) 
k;ar) 


(2.3) 


Proof. Since Xi,X2, ...,Xr are mutually independent, then 

r 

P{Xi<Xi,X 2 <X 2 ,...,Xr<Xr) = WP{Xi < xf) 

i=l 

r 

= l[{l-P{Xi>Xi)). 

i=l 


From fl2.2p . we obtain fl2.3p . □ 

Theorem 2.2. The marginal probability mass function of Xi,i = l,2,...,r is given 
by 


aiMn\l3 + Xj + 1; k; ar) - Mn’{(3 + xp, k; ar) 
M^\(3] k;ar) 


('■)/ 


P{X, = Xi) = 


(2.4) 


Proof. Using 


P{Xi = Xi) = P{Xi > Xi) - P{Xi > Xi) 


and from fl2.2p . we obtain fl2.4p . 


□ 
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2.2 Moments 


The moment generating function is defined by 

Mn\(3; k] 


Mx(t) = 


k] ar) 


(2.5) 


where t = (ti, t 2 , ■■■, U) and X = (Xi, X 2 , Xr). _ 

The £-th moment can be obtained by partial differential of i times the m.g.f in fl2.5p 
with respect to ti and tj = 0 as following 




= iw^x(t)|t,=o= W) 


^_ ^'^r2^r{l — k)n\ 


E 




{x\ -\- X 2 Xi -\-... + Xf + /3)^ 


(n - rk)\M^'^ (d; fc; a^) xi.x2,....x,.....x,=o 
Also, using fl2.4p E^X-) can be written as 

Mn\(5 + Xj + 1; k]ar) - Mn\(3 + x*; /c;^) 


( 2 . 6 ) 


(2-7) 




Xi=0 


Mn\P; k; a) 


( 2 , 8 ) 


Using the method given by Gupta [6], we obtain th factorial moments as follow 
From (12.1 p 




(n — rk)\ 


Xl,X2,...,Xr = 0 


(q(o)^i(q(i)^T..(q;^-i)^" 
(xi +X 2 + ... + x^ + 


the Gth factorial moment can be obtained by partial differential of i times of the 
previous equation with respect to Oj-i as following 


MY\p-,k;a) 




^_j^^r2r(l — fc)n! 


=; E 


{xi)i(ai-iY* ( 01 )""= ...(Qr-i)'' 


(n - rfe)!MM(d; fc; cu) .,..,.,.,(^.....,.0 + ... + + ... + *. + d)’'"- ’ 


( 2 . 9 ) 


where Xi)e = (xi(xi — l)...(xj — I + 1). Since 


i 

{Xi)t = 

j=0 

where s{i,j) are Stirling numbers of hrst kind see |1]. 

Hence, we can obtain relation between moments and factorial moments as follows 


i 


i=o 


( 2 . 10 ) 
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Theorem 2.3. Setting X = {Xi,X 2 ,X^) and t = {ti,t 2 , ■■■, t^) and 
tar = (tiOO) ^ 2 « 2 ) The probability generating function of UGAT distri¬ 

bution is 


Gx(ti,t2, ■■■, tr) 


k] tar) 
k;ar) 


( 2 . 11 ) 


2.3 Bivariate UGAT distribution 

In the case r = 2 in fll.ip for any two random variables Xi,Xj, i ^ j, we have 

2.3.1 Marginal conditional probability distribution 


P{Xi = Xi\Xj = Xj) = 


P ( Xi Xi^ Xj Xj ) 

P{Xi = %) 


(n - 2k)\{xi + Xj)'^’^ ( Mn\f3 + Xi + 1] k] 02 ) - Mn\ld + xp, k] 02 ) 


( 2 . 12 ) 


2.3.2 Marginal expectation 


E{Xi = Xi\Xj = Xj) = XiP{Xi = Xi\Xj = Xj) 

Xi=0 

22(1-A:)^| ^ Xi(ai)^^ 

(n - rk)\ ^ a;^.) 2 fc-n + Xj + 1; fc; 02 ) - + xp fc; 02 )) 

(2.13) 


3 Reliability concepts for UGAT distributions 


3.1 Multivariate reliability function 


Theorem 3.1. Let X = (Ai, X 2 ,..., A^) be a discrete random vector, = {xi,X 2 , ■■■,Xr) 
E representing the lifetimes of r—component system with the multivariate relia¬ 
bility function 


R{x) 


V — 1 

n {at-i)"^^ Mn\/3 + Xi + X 2 + ... + Xr; k-,ar) 

i=l 


k] ar) 


(3.1) 
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Proof. From the definition of the multivariate reliability function, see [12] and fl2.ip 


R{xi,X 2, ■■■,Xr) = PiXi > Xi,X2 > X2: ■■■,Xr > Xr) 

= X! X! ••• X! P{Xi=mi,X 2 = m 2 ,...,Xr = mr) 

mi'>Xi 7712^X2 TTlr^Xr 

{n - rk)\Mi^\f 3 ; k-a,) (^i + ^2 + ... + m. + 

substituting rrii — Xi = it, i = 1, 2,r in the previous equation, we obtain fl3.ip . □ 


Remark 3.1. Let X = (Xi, X 2 ,..., be a discrete random vector, x = (xi, X 2 ,..., x^) 
G Rf. The marginal survival function, see m is defined as 


Ri{ti,X) = P{Xi - Xi > ti\X>x) 

_ R{Xi,X 2, ...,Xi + tj, ...,Xr) 

R{Xi,X 2, ...,Xr) 

SO we obtain the marginal survival function of UGAT distribution 

. Y■^ _ + X 2 + ... + Xj + tj + ... + Xr + /3; k; Ctr) 

(rY~/ ^ ■ 

Mn ( 3^1 + 3^2 + ... + Xr + /?; k] Ctr) 


(3.2) 


Theorem 3.2. R(X) is said to be 


i) Multivariate new better than used (Multivariate new worse than used ) MNBU 
(MNWU) if 


A/W(ft t:;5r)AfP(X, + t + /3; t,ar) < (>)A4’'’(Xr + ft t:;5ftAy>(t + ft ft aft. 

(3.3) 

a) Multivariate new better than used in expectation (Multivariate new worse than 
used in expectation ) MNBUE (MNWUE) if 


r—l 


Mjf\/3-,t,ar) Y1 lla^Mjf\X, + t + (3■,k■,aY < (>) 

tl^t2,...,tr = 0 i=0 


r—l 


MW(Xr + /?;fc;ar) +/3; A:;ar), (3.4) 

tl,t2,...,tr = 0 i=0 

where X = (Xi,X 2 , ...,Xr) be a discrete random vector, Xr = Xi + X 2 + ... + Xr, 
t = ti + ^2 + ... + tr- 
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Proof. When 


MW(Xr + t + /3;A;;a,) < 


+ (5- k-ar)Mlf\t + (5- k-ar) 


k] ar) 


m!{^ (Xr + t + / 3 ; k- ^ M)f> (Xr + /3; k; a,) Mi{> (t + /?; k; a,) 


rip I 


riPi 


Mn\l3; k] ar 


M^\j3]k]ar) Mir\l3;k-,ar) 


ip( R. h.TT, 


+ t + (3] k] ar) ^ 


Mn\l3] k] Ur 


i=l 

r 


ni®-!)”' 


m”(X,+ *:;a,) 


2 = 1 


Mn\j3; k; a 




rj i=l 


. + /?; k] ar) 

k;ar) 


hence, we get 


R{xi +ti,X 2 + t 2 , ...,Xr + tr) < R{Xi, X 2 , ..., Xr)R{ti, ^2, •••, U)- 

From the dehnition of the Multivariate new better used, see [9], then R(X.) is MNBU. 
Similarly, from the dehnition of MNBUE (multivariate new better used in expecta¬ 
tion),see [9], 


y~( R{xi p ti,X2 + t2, •••, Xr + tr) 




R{Xi,X 2 , ..., Xr) 


< ^ R{ti,t2,--,tr) 


? • • • — 0 


we obtain fl3.4p . 


□ 


3.2 Multivariate hazard rate function 

We consider multivariate hazard rate function [12] 

/i(x) = (hi(x), h 2 (x),..., hr{yi)) 

hi(x) = P{Xi = Xi\X > x) 

_ R{pl) ^2; •••) ^i—li I 5 ^i+li ^r) 

R(Xi,X 2, ...,Xr) 

where x = (xi,X 2 , ...^Xr) G R(j_, so we obtain the following theorem. 

Theorem 3.3. The multivariate hazard rate function of UGAT distribution is given 
by 

Mn -f a;2 -t- ... -I- Xi-i -1- Xj -1- 1 -1- Xj+i -1- ... P Xr P Id] k] ar) , , 

hi[yA = l-ai -- (3.5) 

Mk ’{xi P X 2 P ... P Xr P ft] k] ar) 
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Theorem 3.4. Let 

Xr = (Xi + X2 + ... + Xi_i + Xj + 1 + Xj+i + ... + Xr), Xr = (xi + X 2 + ... + x^) an4 

t = + ^2 + ... + 

i/ien the following statement about UGAT distribution holds 
UGAT distribution is MIFR (MDFR) iff 

Mi^\/3 + K + t-k-ar) ^ Mi^\f5 + K-k-ar) ,3 g. 

(/3 + Xr + t; t, Or) {/3 + Xr; k; Or) ' 

Proof, when 


then 


so 


M^^^(/? + X; + t;fc;Or) ^ + X^; fc; Or) 

M^''^(/3 + Xr + t;A:;ar) “ + Xr; A:; Or) ’ 

^ M!f:^ (/5 + X; + t; fc; Or) ^ (/5 + X^; fc; Or) 

1 — OLi - -r~\ - ^ 1 ~ OLi - -r~\ -) 

Mr (/3 + Xr + t; k- Or) Mr (/3 + Xr; fc; Or) 


/l(xi + ti, X2 + ^ 2 , ..., + tr) > h{xi, X2, ..., Xr). 


Hence UGAT distribution is MIFR (multivariate increasing failure rate), see [12]. □ 


3.3 Multivariate mean residual life 

Let X = (Xi, X 2 ,..., Xr) be a discrete random vectors, x = (xi, X 2 ,..., Xr) € R+. 
We consider MMRL(m) [2] 

m(x) = (mi(x),m 2 (x), ...,mr(x)) 


mj(x) = i?(Xj — Xi|X > x) 

00 

= 5^R,(U,X), 

ti=0 

so we get the following theorem 


9 








Theorem 3.5. MMRLm of UGAT distribution is given by 


MMRL{xi) 


OO 


{oii-iY''Mn ^(/3 + Xi + X 2 + + Xj -r (.j -r ••• -r 


ti=Q 


MY\(3 + Xi+ X2 + ... + Xr] k] ttr) 


“1“ tj “1“ ... ^“ Xj.j k] Of.) 


4 Maximum likelihood estimators 


We use the maximum likelihood method to estimate the unknown parameters of 
UGAT distribution. Consider constant values to r, n, we want to estimate the 
parameter at, i = 1,2, ...,r, jS and k. Suppose we have a sample of size N in 
the form ((Xii,X 21 , .■.,Xrl), (Xi 2 ,X 22 , ■■■,Xr 2 ), ..., 

{xin,X 2 n, ■■■,XrN)) from UGAT distribution. Based on observations, the Likelihood 
function of the sample given by 


L{a„/3,k) = l[C 

i=l 


K")(«r)-KLh) 


(xij + X 2 i + ... + XriY^~'^MY\f]; k] Or) 


(4.1) 


{n—rk)\ 


where C = Y\_ 

i=0 

Log-Likelihood function can be written as 


r N N 

logL = logC + ^ log(aj_i)-|-(n-rA:) ^ log(a:ii-|-2:2i-|-...-l-2;ri-|-/3)-A^log(M,('''(/3; k-oir)). (4.2) 

j=li=l i=l 


Computing the hrst partial derivatives of fl4.2p with respect toaj,z = 0,l,...,r — 1, 
{3 and k, we get the normal equations 

N 

N 

- - E{Xi),i = l,2,...,r-l 

ai-i Oii-i 

^ 1 

— (rfc — n) ---h N{rk — n)E{{Xi -|- X 2 -I- ... + Xr + /3)~^) 

^li) “1“ ^2i H“ H” ^ri 

N 

—r ^ ) log(a:ii -|- X 2 i -t... -t Xri + d) + NrE{\og{Xi) -\- X 2 + ... + Xy + (3)) 

i=l 

In the case of r = 3, we study the bivariate UGAT distribution to get the MLEs of 
the parameters cto, ai, 02 and j3. For that we have to solve the above system of four 


dai-- 


■\ogL = 


d , 


dk 


logL = 
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non-linear equations with respect to ao,ai,a 2 and jS. The normal equations are 


^ 1 r 

■ 7 —\ogL 

aao 


dai 


logL 


^ 1 r 


A 

dp 


logL 


N 


S Xij 

- -E{X,), 


i=i 

(4.3) 

ao 

ao 


N 



Xij 

N . ^ 

- —E{X2), 


i=i 

(4.4) 

ai 

Ui 

N 



Xij 

- —E{Xx,l 

Ot2 


i=i 

0.2 

(4.5) 


-{2k - n) V---^-— + N{2k - n)E{{X^ + + P)-^) (4.6) 

^ X2i H“ ~h p 


These equations are not easy to solve, so numerical technique is needed to get the 
MLEs. The approximate conhdence intervals of the parameters based on asymptotic 
distributions of their MLEs are derived. 


5 Data analysis 

In this section, we present the analysis of a real data using the UGAT model and 
compare it with model like multivariate poisson-log normal P /\^ model. The fol¬ 
lowing data represent the study of the relative effectiveness of three different air 
samplers 1,2,3 to detect pathogenic bacteria triplets of a microbiologist obtained 
triplets of bacterial colony counts xi,X 2 , x^ from samplers 1, 2, 3 d the data in each 
of 50 different sterile locations. This data is available in [1] and the data are repre¬ 
sented in the following Table 1. 


Table 1: Bacterial counts by 3 samplers in 50 sterile locations 


Ai 

As 

As 

Ai 

As 

As 

Ai 

As 

As 

Ai 

As 

As 

Ai 

As 

As 

1 

2 

11 

3 

6 

6 

3 

8 

2 

7 

10 

5 

22 

9 

9 

8 

6 

0 

3 

9 

14 

1 

1 

30 

2 

2 

8 

5 

2 

4 

2 

13 

5 

4 

2 

25 

4 

5 

15 

3 

15 

3 

2 

0 

6 

2 

8 

1 

9 

7 

3 

7 

6 

3 

1 

8 

2 

2 

1 

1 

5 

6 

5 

5 

4 

8 

8 

10 

4 

4 

6 

0 

4 

6 

4 

14 

1 

7 

4 

4 

7 

3 

2 

10 

8 

7 

3 

4 

9 

2 

3 

9 

2 

7 

3 

2 

6 

8 

5 

6 

6 

6 

8 

4 

6 

7 

6 

8 

1 

14 

6 

2 

3 

10 

4 

14 

7 

3 

10 

6 

3 

4 

12 

2 

13 

0 

1 

7 

3 

3 

3 

14 

4 

7 

10 

1 

9 

7 

14 

9 

5 

2 

9 

12 

6 

8 

3 

2 

4 

6 


The UGAT model is used to £t this data set. The MLE(s) of unknown parameter(s) 
are Oq = -787, ai = .846, a 2 = .849, fp = 954.707. 


11 













The value of log-likelihood, Akaike information criterion(AIC) and Bayesian infor¬ 
mation criterion (BIC) test statistic two different models are given in Table 2 


Table 2: The MLE(s) parameters, Log-likelihood, AIC and BIC 


The model or hypothesis 

parameters 

-L 

AIC 

BIC 

P/\^ model 

9 

397.8 

813.6 

810.89 

P/\^ with equal fj, 

7 

-402 

818 

815.89 

Equicovariance 

5 

401.1 

812.2 

810.69 

Isotropic 

3 

404 

814 

813.097 

UGAT model 

4 

401.797 

811.594 

810.39 


Table 2 shows that The UAGT model hts certain well-known data sets better than 
P/\^ model. Reducing the number of parameters to four by hxing one of parameters, 
therefore k still provides a better £t than existing model. 
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